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The phase diagram of the XXZ spin-1 quantum magnet on the kagome lattice is studied for all 
cases where the Jz coupling is antiferromagnetic. In the zero magnetic field case, the six previously 
introduced phases, found using various methods, are: the nondegenerate gapped photon phase 
which breaks no space symmetry or spin symmetry; the six-fold degenerate phase with plaquette 
order, which breaks both time reversal symmetry and translational symmetry; the "superfluid" 
(ferromagnetic) phase with an in-plane global U{1) symmetry broken, when J^y < 0; the \/3 x \/3 
order when Jxy > 0; the nematic phase when D < and large; and a phase with resonating dimers 
on each hexagon. We obtain all of these phases and partial information about their quantum phase 
transitions in a single framework by studying condensation of defects in the six-fold plaquette phases. 
The transition between nematic phase and the six-fold degenerate plaquette phase is potentially 
an unconventional second-order critical point. In the case of a nonzero magnetic field along z, 
another ordered phase with translation symmetry broken is opened up in the nematic phase. Due 
to the breaking of time-reversal symmetry by the field, a supersolid phase emerges between the 
six-fold plaquette order and the superfluid phase. This phase diagram might be accessible in nickel 
compounds, BF4 salts, or optical lattices of atoms with three degenerate states on every site. 

PACS numbers: 



I. INTRODUCTION 

The behavior of "frustrated" magnets, in which not 
all interaction energies can be simultaneously minimized, 
is already quite complex when the individual spins are 
treated classically. Models of quantum spins with frus- 
trating interactions are an active subject of current ex- 
perimental and theoretical study. A simple example of 
a frustrated quantum magnet is the standard nearest- 
neighbor Heisenberg antiferromagnet on any lattice with 
closed loops containing an odd number of sites: impor- 
tant examples include the triangular and the kagome lat- 
tices in two dimensions, and the pyrochlore lattice in 
three dimensions. 

For physical magnets with finite values of the spin s, 
there are general approaches such as computing 1/s cor- 
rections to the classical limit s — > 00 and expanding the 
spin algebra from SU{2) to a larger group. Such ap- 
proaches are powerful and predict many interesting or- 
dered phases, but their applicability to real magnets with 
only SU{2) symmetry and small values of the spin (e.g., 
s = 1/2 or s = 1) is uncertain. In recent years, interest 
has shifted to understanding specific examples of finite- 
spin magnets in detail, even though the necessary the- 
oretical methods are less general than either the 1/s or 
large-A^ expansions. Frustrated quantum antiferromag- 
nets with small spin s = 1/2 or s = 1 have been proposed 
to show various exotic behaviors, including gapped or al- 
gebraic spin liquids with gauge-boson-like excitations or 
unconventional second-order phase transitionsii^!^. 

It is often possible to compare such predictions with 
large-scale numerical Monte Carlo studies in cases with 
reduced symmetry (e.g., with SU{2) broken down to 



[/(I)) , but frustrated magnets with full SU{2) symme- 
try are in general accessible only by exact diagonaliza- 
tion, series expansion, or density-matrix renormalization 
group on relatively small systems because of a "sign prob- 
lem" associated with the frustration. The s = 1 model 
on the kagome lattice studied in this paper is motivated 
both by the existence of materials such as BF4 salts^ and 
Ni^+-based materials including Ni3V203^, and by intrin- 
sic interest in the unexpected phases of the model. Our 
goal is to present a single treatment of the two-parameter 
phase diagram of the model that unifies previous studies 
of parts of the phase diagram^'^'- and allows consider- 
ation of the various phase transitions occurring in the 
model. 

Previous theoretical studies on the s — 1 kagome 
lattice antiferromagnet with uniaxial anisotropy ("XXZ 
anisotropy"), with Hamiltonian 

H^Yl -^^^tS^ + D{Stf -f J.,(5f 5J -f SfS^^l (1) 

<ij> 

Here the sum is over nearest-neighbor bonds on the 
kagome lattice. Note that the on-site anisotropy term 
would be forbidden for s = 1/2 and is compatible with 
inversion symmetry, unlike the Dzyaloshinksii-Moriya 
term, also quadratic in spin, that appears if the other ions 
of the crystal break inversion symmetry. For general cou- 
plings, this Hamiltonian breaks the spin rotation symme- 
try SU{2) down to the U{1) subgroup generated by S^, 
and has time-reversal symmetry. We discuss both easy- 
plane and easy-axis limits, and also consider briefly the 
effects of a magnetic field that breaks time-reversal but 
preserves the t/(l). Section II reviews previous theoreti- 
cal work on the zero-temperature physics of this Hamilto- 
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nian, which for different values of the couplings has found 
a gapped phase with a massive photon-like excitation^, 
a critical line separating plaquette-ordered phases^, and 
an Ising-type spin ncmatic^. Section III presents the 
field-theory description of the plaquette-ordered phases 
in terms of dual height variables. From Section IV to 
section VII, we study the transitions between the six- 
fold degenerate phase and other phases, we will see that 
all the other phases can be interpreted as the conden- 
sates of different kinds of defects in the six- fold degener- 
ate plaquette phases. In section VIII, the situation under 
longitudinal magnetic field (along z) is studied, several 
new phases are found. Section IX is devoted to the point 
with spin-S'[/(2) symmetry, and section X is about other 
transitions in phase diagram Fig. O 



II. EXPERIMENTAL SYSTEMS AND 
PREVIOUS STUDIES 

So far two types of kagome spin-1 materials have 
been found. The first type is BF4 salts^, the second is 
Ni^+ based material Ni3V2 085. Also, the kagome lat- 
tice has been constructed with laser beams-, an effec- 
tive spin model can also be realized in cold atom sys- 
tem trapped in optical lattice, but there the existence of 
biquadratic interactions comparable in strength to the 
standard Heisenberg interaction makes the phase dia- 
gram even more complicated^^. 

In general the model we are going to discuss is de- 
scribed by equation ([1]) . This Hamiltonian is the simplest 
example which can potentially realize all the physics dis- 
cussed in this work, but our formalism is supposed to be 
more general, and independent of the details of the model 
on the lattice scale. This is the simplest spin model 
which is invariant under time reversal transformation. 
Three coefficients , J^y and D are used to parameterize 
this model. If all the coefficients are positive, this model 
can be realized in magnetic solids like those given above; 
when Jz,D > 0, and Jxy < 0, this model could possibly 
be realized in cold atom systems with pseudospin degrees 
of freedom on each site. For instance, suppose on every 
site there are three orbital levels (the three orbital levels 
can be the degenerate p-level states, as discussed in sev- 
eral previous papersii), the orbital degrees of freedoms 
can be viewed as spin-1 pseudospin, with natural XXZ 
symmetry. The antiferomagnetic coupling and D can 
be generated by the on-site s-wave scattering and off- 
site dipole interactions^^' The J^y coupling is resulted 
from the superexchange, which should be ferromagnetic 
due to the bosonic nature of the system. Therefore in 
the following discussions, both positive Jxy and negative 
Jxy cases will be discussed. 

In solids, the spin SU{2) symmetry can be broken by 
spin-orbit coupling and the layered nature of the mate- 
rial, or by an external magnetic field; in the cold atom 
pseudospin system, the SU{2) symmetry is missing at 
the very beginning, as the orbital level pseudospin sys- 



tem has natural uniaxial anisotropy. 

Several previous papers have studied the kagome spin-1 
system^i^iii, at different parameter regimes of this par- 
ticular model ll]). There are five phases that are already 
known. 

Superfluid phase: When Jxy < and \Jxy \ S> D,Jz, 
in this case Jxy is the dominant term in the Hamilto- 
nian ID): the expected phase is a superfluid phase that 
breaks the global U{1) symmetry. In the spin language 
this phase is a ferromagnet in XY plane. Here the term 
superfluid phase is used since the broken symmetry of 
this phase is the same as the superfluid phase. 

\/3 X Vs phase: When Jxy > 0, and becomes the 
dominant term in the Hamiltonian, the phase is not ob- 
vious at first glance. When D — Jz — 0, and the spin 
S +00, the system is at the classical XY limit. It has 
been shown that the ground state of this classical XY 
model has a large discrete degeneracy, in addition to a 
U{1) that rotates all the spins: the zero-temperature en- 
tropy associated with this degeneracy is proportional to 
the size of the system. The ground state configurations 
satisfy the requirement that every triangle has zero net 
spin. If one spin is fixed, the whole ground state configu- 
rations can be one to one mapped to the classical ground 
states of the three-color modeli^. Three-color model is 
defined as follows: on the honeycomb lattice, each link 
is filled by one of the three colors, green, red and blue, 
and the whole lattice is colored in such a way that every 
site joins links of all three colors . The classical parti- 
tion function is defined as the equal weight summation 
of all the 3-color configurations. This partition function 
and entropy have been calculated exactly by Baxter—. 
It has also been shown that the classical model can be 
mapped to a critical 2-component height model (similar 
to our model) ^^'^^ , the low energy field theory of this 
model is a c = 2 conformal field theory with SU{3)k=i 
symmetryiSii^. 

The large degeneracy of the classical model is not uni- 
versal, and it can be easily lifted by the second and third 
nearest neighbor interaction J2 and J3. When J2 > J3, 
the q = state (Fig. [2]) is stabiHzed; while if J3 > J2, 
the \/3 X \/3 state (Fig. [T]) is stabiHzed^a. The large 
3-color degeneracy is also lifted by 1/S expansion, and 
some ordered pattern is picked out from all the classi- 
cal degenerate ground states, this effect is usually called 
"order from disorder". At the isotropic case (J^ — Jxy, 
D — 0), it was proved that after 1/S expansion both 
coplanar q — state and the \/3 x \/3 state are stable 
i.e. they are both local minima in all the ground states, 
the spin wave modes around these two minima do not 
destabilize the order. Latter on, more detailed studies 
suggest that the global minimum state is the -s/S x -s/S 
order—, as depicted in Fig. [T] Although the 1/S expan- 
sion is carried out at the isotropic point, the coplanar 
-y/S X phase is expected to extend to the limit when 
Jxy is dominant. 

Gapped photon phase; When \Jxy\ <C Jz < a 
gapped phase without any symmetry breaking has been 



FIG. 1: the ^ X ^/3 order. 
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FIG. 2: the (7 = order. 



found^. The low energy excitation with the smallest gap 
is a loop excitation with the same polarization and gauge 
symmetry as a photon: the effective theory can be de- 
scribed as a one-component massive compact gauge field. 

Plaquette phase: When \Jxy\ ^ D and jj^ — 13 1, 
< D < Jz, a gapped phase with a six- fold degenerate 
ground state has been found^. The six- fold degenerate 
ground state has plaquette order: spins resonate around 
a subset of the hexagons in the kagome lattice. In this 
parameter regime, the classical part of this model can be 
written as 

^ - E y(E stf + - J^)(sn'. (2) 

A i=l i 

When < D < Jz, the classical ground states are 
all the configurations with every triangle occupied by 
— (1,-1,0). Again the classical ground states can 
be mapped onto the classical 3-color model^, although 
the 3-color states correspond to instead of spins in the 
XY plane (Fig. E]). 

If small Jxy is turned on (either J^y > or J^y < 0), 
the large degeneracy of 3-color ground states is lifted, 



FIG. 3: The classical ground states of ([T}, when J^y — 0. Ev- 
ery triangle has configuration (1,-1,0), which can be mapped 
to the 3 color model on the dual honeycomb lattice. 

and the effective Hamiltonian which operates on the low 
energy Hilbert space is 

H^J2 -tiSlS^SlS^SlS^ + H.c). (3) 
o 

1 to 6 are the sites of each hexagon on the kagome lattice. 
The flippable hexagons have four kinds of configurations, 
they are (1,0,1,0,1,0) (denoted as Ai), (0,1,0,1,0,1) 
(denoted as ^2), (—1, 0, —1, 0,-1, 0) (denoted as Bi) and 
(0, — 1, 0, — 1, 0, — 1) (denoted as 82)- The ring exchange 
term ([3]) can flip Ai to A2 (and vice versa) Fig. IH also 
can flip Bi to B2 (and vice versa). Two compact U{1) 
gauge fields were introduced to describe this system, and 
due to the monopole proliferation, the system is gener- 
ally gapped, with crystalline order. The particular order 
which happens here is the plaquette order, which breaks 
both translational and time-reversal symmetries. The 
simplest way to view this state is that, since the ring ex- 
change term ([3]) can flip either Ai to A2 configurations, 
or flip Bi to B2 configurations, the configurations with 
the largest number of flippable hexagons are favored in 
order to benefit from this ring exchange term. Notice 
that the hexagons form a triangular lattice with three 
sublattices, then one out of the three sublattices of the 
triangular lattice can be resonated. Also one can choose 
either to resonate between Ai and A2 configurations or 
to resonate between Bi and B2 configurations (these two 
resonance cannot both happen at the same state). There- 
fore there are in total 3x2 = 6 degenerate ground states. 

The simple picture of the ground state will be fur- 
ther justified in the next section, by studying the dual 
quantum height model. The classical height model was 
introduced to study the classical 3-color model, and since 
there are two components of free boson height fields in 
the continuum limit, it is believed that the low energy 
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FIG. 4: The effect of the ring exchange term It can flip 
Ai [Bi) to A2 (S2), and vice versa. 



field theory should be c = 2 CFTJ2. We will see that the 
quantum effect is relevant at the classical 3-color critical 
point, a gap is opened due to the vertex operators of the 
height fields. 

Recently a mean-field treatment of a similar model has 
been studiecU^. The plaquette phase we obtained is sim- 
ilar but not entirely identical to the "Plaquette ordered 
phase" in this recent work, which is identified as the fully 
packed string crystal. The main difference between the 
two approaches is that, the monopole effect of compact 
gauge theory has been taken into account in our work 
from the very beginning. The monopole effect is sup- 
posed to be very relevant at the z = 1 Gaussian fixed 
point of gauge theory, and dominate the physics close to 
the Gaussian fixed point. The nonlocal monopole effects 
can be described by a local field theory in the dual for- 
malism, and the ordered pattern is predicted in this dual 
local field theory. 

Nematic phase: When D < and \ Jxy\,Jz <C \D\, 
a nematic phase with nonzero expectation value of (S'^)^ 
has been found*. In this case, because D is negative and 
large, the system favors S"^ = ±1 on every site. Since 
the = state costs too much energy, every site can 
be viewed as an Ising spin, and this model is effectively 
equivalent to a spin-1/2 model. Since {S~)'^ flips = 1 
to 5^ = — I state, it plays the same role as = — ia"^ 
on the effective Ising spin. Therefore the superfluid phase 
of this spin-1/2 system is actually the nematic phase of 
the original spin-1 model. 

The rough phase diagram is shown in Fig. [Sj The 
goal of the current work is to understand all the phases 
we know from the excitations of the six-fold degenerate 
phase. Basically all the phases can be interpreted as the 
condensates of various defects which violate the (1, —1, 0) 
constraint in the plaquette phase above. Since the low 
energy Hilbert space is a constrained one, to create one 
single defect cannot be realized from local moves of the 
ground state configurations, instead, global change of all 
the spins is involved. This implies that one defect in 
this phase not only carries the global U{1) charge, but 
also carries the gauge charge, with the gauge symme- 
try emerged at the low energy constrained Hilbert space. 
Therefore the condensate of defects is also the Higgs 
phase of the compact gauge fields. 




Disorder 



FIG. 5: the phase diagram with zero magnetic field. When 
D > Jz and \Jxy\ is small, the system is in the nondegen- 
erate gapped photon phase; when < D < Jz, and with 
small \Jxy\, the ground state breaks both translational and 
time reversal symmetry, resulting in the six-fold degenerate 
plaquette order; when J^y is negative and large, basically the 
system is in superfluid order which spontaneously breaks the 
global spin conservation; when Ja;j, is positive and large, 
X y/3 order is supposed to be favored; and when D is neg- 
ative and large, the system has nematic order, with nonzero 
expectation of (S^)^. 



III. GAUGE THEORY OF THE 
PLAQUETTE-ORDERED PHASE 

When = and < D < J^, the set of degenerate 
ground states can be mapped exactly! to those of the 
3-color modeli^. Every triangle on the kagome lattice 
has configuration (1,-1,0) on this classical critical 
line. The z-component spin configuration on the kagome 
lattice can be viewed as two-component dimer configu- 
rations on the dual honeycomb lattice, with repulsive in- 
teraction between two flavors of dimers (every link of the 
honeycomb lattice can only be occupied by one dimer). 
It is well-known that the one component quantum dimer 
model can be mapped to compact gauge theory^^, there- 
fore it is natural to describe this spin-1 system as two 
compact U{1) gauge flelds, since we can interpret the 
3-color constraint as two independent U{1) constraints: 
every site on the honeycomb lattice connects to exactly 
one — 1 dimer and one S"^ = — 1 dimer. We may 
map the three values of to three configurations of a 
two-component electric field: 



5f = -1 ^ iE„E,) ^ i=(-i, -^). (4) 
Next, note that a 2D unit vector can be assigned 
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parallel or antiparallel to each link i of the honeycomb 
lattice (dual lattice of the kagome lattice) so that vertices 
of sublattice A of the honeycomb have three incoming 
bonds, while those of sublattice B have three outgoing 
bonds. Now define two-component vector fields on bonds: 
Ea = Eaiii. The three color constraint is now equivalent 
to the Gauss's law constraint 



consider the defects, which are also the gauge charges. 
When D is much smaller than Jz — D, the excitation 
with the smallest gap is to flip one site with 5^ = to 
1 (or -1), this process actually creates a pair of (1,1,-1) 
(or (—1, —1, 1)) defects. Let us denote the density of the 
(1,1,-1) configuration defect as pi, and denote the den- 
sity of the (1,-1,-1) defect as p2, then from the defini- 
tion of electric field we can obtain the following relations 



(5) 



Also we can generalize the configuration of the E vec- 
tor to a 2d triangular lattice. The lattice is formed by 
basis bi = (\/3/2, 1/2) and h = (0, 1), 



E^nbi+ mh - (l/(2\/3), 1/2) 



(6) 



If we add the following interaction to the Hamiltonian, 
the E fields only take three smallest vectors as Q: 
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E-' 



(7) 



Thus the low energy configurations of electric fields can 
be one-to-one mapped to the low energy configurations of 
spins; the spin formalism and the electric field formalism 
are equivalent. 

The perturbation theory of Jxy generates a ring ex- 
change (III). The ring-exchange term breaks the Z3 
symmetry. Define conjugate operators on each bond 
(Ai^i, A2^i) with commutation relations 



(8) 
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P2) 



(11) 



The charges can be effectively viewed as matter fields 
defined on the sites of the honeycomb lattice, and the 
gauge fields A and E are fields defined on the links of the 
honeycomb lattice. 

For the convenience of later calculations, we need to 
define a new set of variables as follows 



1 



ei = -El + E2 

v3 

r 1 7 
ai = --M + -A2 

^ A It 

a2 = ~—Ai--A2. 



(12) 



Also, one can check that Ca and (Jq, are still conjugate 
variables : 



Then operator Ta.i = eyi\>{iAa.j) acts as a raising oper- 
ator: it increases the quantum number Ea.j by 1. This 
enables a compact representation of the ring-exchange 
terms proportional to t: on bond j, exp(iy4.a jZi^^') will 
raise SI = to SI 



larly, if 



1 if = (-V3/2,l/2). Simi- 
= (-\/3/2,-l/2) then exp(i^Qj/i^^) takes 

the 



= to 5| = —1. Define vector A^^j = A^j-nj, 
ring exchange term around each hexagon becomes 



-tcos(V X AJ^^^) 



(9) 



(13) 



If the definition for Cq, and a_a is plugged in pip , one 
can see that is the electric field corresponding to the 
charge pa, in the sense that 



V • = -pa- 



(14) 



When D is smaller than but close to Jz, \Jz — D\ <^ D, 
the lowest energy excitation is (0, 0, 0), and we denote its 
density as po- It carries gauge charge of gauge field Ei 



Here as usual in gauge theories of lattice spin models, 
the meaning of V x ^ is that one takes the lattice circula- 
tion around a plaquette: for Vj an clockwise assignment 
of unit vectors along the links around a hexagon 



A, 



(10) 



Po 



V3 



{V -ei+V ■ 62) 



(15) 



Since the electric fields are subject to the constraint 
(ini), it is convenient to define height fields h on the dual 
triangular lattice. 



If no defect is present, i.e. the Gauss's law constraint 
is strictly imposed, the theory is described by two com- 
pact U{1) gauge fields without matter fields. Now let us 



Ea^ {z X V)/la, V X Aa= TTha- 



(16) 



TTfia and ha are a pair of conjugate variables. The value 
of ha is also defined on a triangular lattice configuration 
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space, in order to satisfy definition (|T6l) . ha is defined in 
the following way 



V3 1 

{hi,h2)a = (— (m + n),-(TO-n)) + {qi,q2)a, (17) 

m and n are both integers. Here a — A, B, C, denoting 
the three sublattices on the triangular lattice (dual lattice 
of the honeycomb lattice), qa are three vectors, taking 
different values on three sublattices 



= (- 



1 



^2\/3' 6 
^^^(O i) 



(18) 



The two-component height variables are the same 
as those introduced in the classical 3-color model (cf. 
Kondev and Henlej^). Since m and n are both inte- 
gers, the vertex operators should enter the effective low 
energy theory. We will see later that, due to quantum 
effect, these vertex operators become very relevant and 
drive the system away from the classical criticality, re- 
sulting in a six- fold degenerate plaquette ordered phase. 
These vertex operators read 



1 1 

Hv = -aC0s[27r(— =/lia + /l2a + -JElla + q2a)] 
a=A 

- aCOs[27r(— =/lia - h^a + —f^lla - g2a)]-(19) 
a=A V3 V3 

For later convenience, we define a new height fields 
and its conjugate variable Tr^a as 



/l2 
/»2 



7r2 



V3 

2 
%/3 



2 ''^hl + '^''^hl 



(20) 



One can check the commutators and see that (jja and 
are conjugate variables, and based on the definition 
they are exactly the height fields corresponding to 
e and a. 



{z X V)0Q, V X da 



(21) 



The vortex of (pa is the charge field pa- 

Now in terms of the new height fields, the vertex op- 
erators read 



2tt 

Hy = -a[cos(27r(/)i^) -I- cos(27r(/)iB + — ) 

o 

-hcos(27r0ic + — ) 
47r 

-I- cos(27r(/)2A + ^) + cos(27r(?!)2s) 

27r,, 

-f cos(27r02C + "3")J- 



(22) 



These vertex operators have oscillating signs on the 
triangular lattice, then in the low energy theory the rel- 
evant terms should be higher orders of vertex operators 
which do not contain oscillating signs on the lattice: 



= —v[cos{6'Kipi) + cos(67riy92)] 
-ui [cos(27ri^i + 4n(p2) + cos{4TTipi + 2inp2)] 



-V2 cos(27riy9i — 27r(^2) + 



(23) 



ipa is the coarse-grained mode of (j^a- As the vertex 
operators correspond to the creation and annihilation of 
gauge fiuxes, the total gauge flux is conserved by mod 3 
in the low energy continuum limit. In equation (|23p . v 
is supposed to be positive, but vi and V2 are supposed 
to be negative, because when we subtract (j>2 from 0i 
from (im, it gains angle 27r/3, which generates a factor 
— 1/2 before the cosine term in ((23)l . Sine functions of ipa 
are excluded by symmetries of the system. For instance, 
sin(27r(93i — 1^2)) is excluded by time reversal symmetry. 

After coarse-graining the system, the action in terms 
of ipa can be written as 

2 

^ = Zl^^^"^")^ + P2(Vv?„)^ + Hy + 7V(^iV(^2 (24) 



The 7 term in ()24|) is a flavor mixing term between tpi 
and (p2 , and therefore the two flavors of height fields do 
not only couple to each other through the vertex opera- 
tors, but also through one of the kinetic terms. 

In 2-|-ld, the potential operators with cosine functions 
are generally very relevant at the fixed point described 
by the Gaussian part of the action as long as the 

term {p2 in ([24]) ) is present. Vertex operators are 
responsible for the gapped crystalline phases of quantum 
dimer models, both on the square lattice^ii^^ and the 
honeycomb lattice^. In the current work, the vertex 
operators are also responsible for the crystalline phases. 
First of all, let us tune vi and V2 to zero, and minimize 
V terms in (|23p . Each tpa has three minima , 1/3, 
2/3. Therefore there are in total 9 different combinations. 
However, negative vi and V2 terms will raise the energy of 
all the minima with (pi — (f2, and hence we end up with 
9—3 = 6 minima. This result is actually quite general, for 
a large parameter regime, there are always 6 minima of 
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the vertex potential in ((23)) . Because the vertex operators 
in (|23p is invariant under transformation ipi —^ ipi + 1/3, 
(y52 ^ + 1/3, and also invariant under transformation 
ipi ^ ip2, all six minima can be obtained from performing 
transformations on one single minimum. 

Now we can write down the plaquette order parame- 
ter in terms of the field theory variables (fa- The order 
parameter we are searching for, in the lattice model, is 

a=A O 

(25) 

In the above equation, a represents different sublat- 
tices, and = 1, is = 2, and ic = 3. The low energy 
representation of this order parameter can be deduced 
from symmetry argument. The most obvious transfor- 
mations for this order parameter are translational (T) 
and time reversal (TR) transformations. 

T : e*27r/3p^ TR : P ^ -P. (26) 

If rotated around hexagons at sublattice A by angle 
27r/3 (R27r/3)i the order parameter is invariant; under 
space inversion (SI) f—^ — r and reflection (Px) along x ( 
y —f —y ) centered at sublattice A, the order parameter 
becomes its complex conjugate 

SI,Px:P^P*. (27) 

Under the transformations discussed above, ipi trans- 
forms as follows: 

T : ipi ^ ipi + ^, 
TR : <y9i ^ ip2, 

SI: tpa -fa, 
Px : V'q ^ -fa, 

R27r/3 : fi fi- (28) 

Summarizing all the transformations above, the field 
theory representation of the plaquette order parameter 
P is 

p ^ (^^^^^2 _ ^i2^vi y (29) 

We can plug in the six minima of the vertex opera- 
tor to and it gives us 6 different values. All 
the six expectation values can be obtained by transfor- 
mation (Pn) = exp(i7rn/3)(Po), with n = 1 to 6. This 
implies that the system is in a plaquette order with six 
fold degeneracy. 

When both vi and V2 are positive, the vertex operators 
gives three degenerate minima: fi^f2 — 0, or ±1/3. 
These three degenerate ground states do not break time 
reversal symmetry, but it breaks translational symmetry. 
The particular order in this case is another type of pla- 
quette order with (—1, 1, —1, 1, —1, 1) hexagons resonat- 
ing on one of the three sublattices. 

The field theory description is only valid when the the- 
ory is close to a critical point, i.e. the correlation length 



is either infinite, or finite but much longer than the mi- 
croscopic lattice constant. Thus the prediction of pla- 
quette order is only rigorous close to the classical critical 
point Jxy = 0. But the phase is expected to extend over 
a finite region in the phase diagram, until a transition 
into either a disordered phase or one of the other ordered 
phases derived in the following sections. 

In this section we started with mapping the classi- 
cal ground states of the model onto the classical 3-color 
model configurations, as in this way we respected the Z3 
symmetry of the classical ground state, which is broken 
by the quantum perturbation. As mentioned before, we 
can also view the low energy physics of this system as 
two components of quantum dimer model, with repulsive 
interaction between two flavors of dimers. From this ap- 
proach the same low energy action as (p4|) can be derived. 
Single component of quantum dimer model generates the 
kinetic terms and the vertex operator cos(67r(/3i) in (j24p 
at low energy, as discussed in referenceSS; the repulsive 
interaction between the two flavors of dimers will gener- 
ate the term ^VipiVf2 and the mixture vertex operators 
cos(27r((^i — 1P2)) and cos(27r(pi -|- 47r(/32), et,al. 



IV. TRANSITION TO THE FEATURELESS 
GAPPED PHOTON PHASE 

When D > Jz, the classical ground state is — 
on every site, and the low energy excitations are 
(-1-1, — 1, -1-1, — 1, ....) loops. This phase has a single 
ground state and gapped photon excitations-, without 
any symmetry breaking. In this section we are going to 
study the phase transition between the six-fold state and 
the gapped photon phase. 

If we start with the six-fold degenerate phase, the tran- 
sition can be viewed as condensation of (0, 0, 0) defects. 
The gap for (0, 0, 0) defect keeps decreasing as the tran- 
sition to the gapped photon phase is approached. But 
the phase boundary between the plaquette phase and the 
nondegenerate phase is not exactly at D ^ (Fig. O, 
this is due to the fact that at second order perturba- 
tion of Jxy/ Jz, an additional nearest neighbor diagonal 
interaction is generated. (0, 0, 0) triangles are more favor- 
able than (1,-1,0) triangles to this diagonal term gener- 
ated, therefore the second order perturbation effectively 
increases D by ^ J^y/Jz- 

The defect (0,0,0) carries charges of both aifj, and 
a2fj., and defects at different sublattices of the honey- 
comb lattice carry opposite gauge charges. If we denote 
the (0, 0, 0) defect at sublattice A as tpA and the (0, 0, 0) 
defect at sublattice B as ipB, the effective Lagrangian 
describing the system close to the transition is 

L = -t\{d^ - iai^ - ia2^)?/'Ap - t\{df_, + iai^ + M2^)i/'sP 

(30) 

Notice that the (0,0,0) defect carries zero global U{\) 
charge (a (0,0,0) defect does not carry any S^), and 



therefore one 4' A particle and one ■i/'s particles can be an- 
nihilated together, so the term g{'ipAtpB + H.c.) is allowed 
in the interaction. After the condensation of ipA and ipB, 
the gauge field = oi^ + will be gapped out along 
with the phase mode 9a — Ob, and the mode 9a + 9b 
will be gapped out by the interaction g{ipAipB + H.c.) 
{9a and 6b are phase angles of V'yi and tps respectively). 
Therefore in the condensate there is no gapless excita- 
tion, which is consistent with the gapped photon phase. 

To further justify this picture, let us first take a 
Landau-Ginzburg tour to study this transition. Let us 
define complex field <& to describe the low energy mode 
of the plaquette order parameter P. The LG action for 
this transition is as follows 

L = - r|$p + + + H.c). (31) 

Without the g term, the theory describes an 3D XY tran- 
sition. The g term turns on an Zg anisotropy at this 
critical point. In the ordered state of this anisotropy 
is a relevant perturbation and will lead to a six-fold de- 
generacy. At the 3Z? XY critical point, Zg anisotropy 
is irrelevant, thus the Landau-Ginzburg theory predicts 
that the transition between the six-fold states and the 
featureless gapped photon phase is a 3D XY transition. 

The 3D XY transition is driven by the vortices of $, 
and after the condensation of the vortices, the vortex core 
state grows and becomes the macroscopic order. It has 
been shown before that the vortex core of the height field 
of the quantum dimer model is an unpaired spin, which 
implies that the condensate of vortices breaks the spin 
SU{2) symmetry spontaneously (for instance, the Neel 
state). In our case, the vortex configuration of $ (includ- 
ing the core) has been depicted in Fig. [HI Around every 
vortex core, there are 6 domains separated by domain 
walls, each domain is one state out of the six-fold de- 
generate plaquette ordered states. In the ordered phase, 
the vortices are linearly confined due to the pinning po- 
tential <1>^ -I- H.c.^ because the domain walls would cost 
energy proportional to their length. At the critical point 
since the pinning potential is irrelevant, the vortices are 
deconfined. 

In Fig. [HI one can see that the vortex core is actually 
a (0,0,0) triangle, which is the lowest energy defect when 
D ~ Jz. If the height field representation of $ ([M]) is 
taken, one can see that the vortex of $ is a bound state 
of one vortex of tpi and one vortex of • Thus a vortex of 
$ carries one gauge charge of ai^ and one gauge charge of 
a2^, i.e. this vortex carries the same gauge charge as the 
(0,0,0) defect. Therefore indeed the transition between 
the six-fold plaquette state and the featureless photon 
phase is driven by the (0, 0, 0) defects. 

The dual field theory of pTjl would describe the vor- 
tex condensation directly. After the standard superfluid- 
gauge field duality in 2+ Id, the dual theory reads 

L = -t|(a^-zA^)^|2 + ... (32) 

Herein ^ is the vortex creation operator, and the Zg 
anisotropy term in (I3ip becomes the monopole processes 
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FIG. 6: The vortex and vortex core configuration in the 
six- fold plaquette order. Every vortex is surrounded by six 
phase domains and 6 domain walls, the vortex core is ex- 
actly a (0,0,0) defect. In this figure, B-|- denotes resonat- 
ing (1,0,1,0,1,0) plaquette on sublattice B. The six do- 
mains around this vortex core are (count counterclockwise): 
with $ ~ $0, and 1^ = 0o + (1/3,2/3); C- with 
$ ~ exp(i7r/3)$o and <^ = (^0 + (1/3,0), B+ with $ ~ 
exp(z27r/3)$o and = i^o + (2/3, 0); A— with $ ~ exp(i7r)<l>o 
and 0^00 + (2/3,1/3); C-|- with $ ~ exp(j47r/3)$o 
and (p — (fo + (0,1/3), B— with "1> ~ exp(i57r/3)$o and 
<^ = (^0 + (0, 2/3). The arrows show the circulation of the 
phase of $. Also, this vortex is a bound state of one vortex of 
ifi and one vortex of (fi2. The spins are defined on the links of 
the honeycomb lattice, which are the sites of the kagome lat- 
tice, links with -I- and - are occupied by spin states S'^ — +1 
and —1. The red dashed lines denote the domain walls. 

which annihilate and create the fluxes of gauge field A^. 
By comparing equation ([5^ and equation (PO]) . we can 
see that A^^ = aifj. + a2fj., and i/j = V'A = V'b- Please note 
that because ipA and V'b can annihilate together, there 
is actually only one flavor of defect, and tpB = i'A- In 
the ordered phase the gauge field is gapped out by 
monopole proliferation, and tp is confined; in the nonde- 
generate photon phase the gauge field is gapped out with 
Tp through the Higgs mechanism. The gauge field is only 
gapless at the critical point. 

If we plug in the height field representation of $ in 
equation ([M)) into the LG action (PT|) . it reproduces the 
height field action in equation , thus the phase transi- 
tion between the gapped photon phase and the plaquette 
phase can also be studied in the dual height model. We 
define new height fields ip± ~ (tpi ± ip2)/2, and they sat- 
isfy the following relation 

(ei ±e2)/2 = (z X V)(/?±. (33) 
Now the height field Lagrangian reads 

L = tidrip+f + P2{Vip+f + t'{drip-f + P2(V(^-)' 



—V cos(67r(^+) cos(67r(p_) — vi cos{6TTip^) cos{2TTip^) 

—V2 cos(47r(p_).(34) 

One (0,0,0) defect carries one unit gauge charge of 
(ei + 62)72, thus it is one unit vortex of ip+, and the 
condensation of (0, 0, 0) drives (p+ into disordered phase. 
In the condensate, (/3-|_ is disordered and the expectation 
value of cos(27r(^_|-) is zero. Thus the plaquette order 
parameter 

P ^ e^^^V2 _ gi27r.pi isin(27rv3_)exp(i27r(/j+) (35) 

takes zero expectation value: the plaquette order disap- 
pears. Since ip- does not transform under translation or 
rotation by 27r/3 transformations, any crystalline pattern 
which breaks these symmetries cannot exist. 

When field (p^ is disordered, the ordered pattern and 
symmetry of the ground state can be studied from the ef- 
fective action for height field ip^ , which remains gapped 
and ordered. Thus the order of the condensate is de- 
termined by the series of vertex operators of if-, since 
the leading vertex operator is —V2 cos(47rc/5_), for a large 
range of parameters, the minima are at ip^ = ±1/4. 
However, let us imagine writing down a physical order 
parameter which only involves ip- — {ipi — ^2)/2, since 
any physical order parameter should be invariant under 
transformation ipa — *■ (pa + i, this order parameter should 
be invariant under (p- — > ip^ + l/2, thus the ground states 
(/?_ = ±1/4 are physically equivalent to each other, and 
the ground state is nondegenerate, which is again consis- 
tent with the gapped photon phase. 

At the transition, since only ip^ is ordered, we can plug 
in the minimum of ip^ into p4p . and obtain an effective 
action for ip^ . Notice that both v and vi vertex operators 
vanish after plugging in the minima ip- = 1/4. The 
leading operator that survives is cos(127r(y5+), which is a 
Zg anisotropy. The height field theory which describes 
this transition is 

L = {d^,Lp+)^ ~ acos{12TT(p+). (36) 

This action describes an XY transition as the Zg 
anisotropy term is irrelevant at the XY critical point. 
Thus we conclude that the transition between the six- 
fold state and the featureless photon phase is driven by 
the condensation of (0, 0, 0) defect, and the critical point 
belongs to the 3D XY universality class. 

V. TRANSITION TO THE SUPERFLUID STATE 

When J^y is negative and large, the system is in the 
superfluid phase (ferromagnetic phase in spin XY plane), 
with nonzero expectation value of (5'^). In this section 
we are going to study the transition between the six- 
fold degenerate plaquette phase and the superfluid phase. 
Let us first focus on the region where D \Jz ~ D\; in 
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FIG. 7: The vortex and vortex core configuration of ipi. Ev- 
ery vortex is surrounded by 3 phase domains and 3 domain 
walls, the vortex core is exactly a (1, 1, —1) defect. The three 
domains around this vortex core are (count counterclockwise) : 
B+ with 0= ipo + (2/3, 0); B- with ip ^ (po + (0, 2/3), A+ 
with (p — ifo + (1/3,0). The arrows show the circulation of 
'Pi- 



this parameter regime, the defects with the lowest gap 
are (1, 1, —1) and (1, —1, —1) triangles. It was shown in 
section III that these two defects are vortices of ipi and 
ip2 respectively. As an example, a vortex of height field 
tpi is shown in Fig. [71 one can see that the core of this 
vortex is a (1,1,-1) defect. When the vortices of the 
height fields condense, which means the height fields are 
disordered, the system enters a superfluid phase. When 
D and J^y are small, the phase transition occurs when 
the hopping energy of the defects is comparable with the 
gap, the phase boundary is roughly D ^ J^y, as shown 
in the phase diagram Fig. O 

Defect (1,1,-1) can stay at two sublattices of the hon- 
eycomb lattice, let us denote defect (1,1,-1) at sublat- 
tice A as V'lA, denote defect (1, 1, —1) at sublattice B as 
i'lB', denote defect (1,-1,-1) at sublattice A as '(/'2A, de- 
fect (1, —1, —1) at sublattice B as '!/'2B- Herein 4 flavors of 
defects are defined because these defects have indepen- 
dent conservation laws instead of just one global U{1) 
conservation law in the original Hamiltonian. If we want 
to hop one (1,1,-1) defect from sublattice A to sublat- 
tice B, global spin configurations within the low energy 
subspace should be changed; this means that any local 
operator cannot hop (1,1,-1) defect from sublattice A to 
B, i.e. defects at sublattice A and B are separately con- 
served. This situation is similar to the doped quantum 
dimer model on square lattice^^, in that case the doped 
holes can also only hop in one sublattice due to the gauge 
symmetry of the dimer model. The gauge symmetry of 
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the dimer model is due to the dimer constraint imposed 
automaticaUy. 

According to equation dU]), ipiA (V'ib) carries charge 
+1 (-1) of gauge field ai^, and V'2A (V'2b) carries charge 
+1 (-1) of gauge field 02^. Now the effective Lagrangian 
describing the system is 

L = -t\{df, - iaif,)'ipiA\'^ - t\(df, + Wi^)i/'ibP 

-t\{d^ - ia2p)V'2yiP - t\{d^ + ia2p)ii2B? 

+ • (37) 

The ellipses include the monopoles of gauge fields as 
well as the interaction terms between different matter 
fields. The interaction has to be consistent with all the 
internal symmetries of the system, which is U(l)giobai x 
U{l)gauge X U{l)gauge- The U {1) gaugc Symmetries cor- 
respond to the two flavors of gauge fields, and the U{1) 
global symmetry corresponds to the conservation of S^. 
The regular terms like — r|'(/'Mp + O(|'0|^) are all allowed, 
besides these terms, another term should in principle ex- 
ist, which is —g[4'iA'4'iB'4'2A'4'2B + H.c ). Four different 
flavors of particles can be created and annihilated to- 
gether, without any global reconfigurations. 

The superfluid phase can be viewed as the condensate 
of 4 flavors of matter fields. Let us denote V'm ~ exp(6'ia), 
the action can be written as 

L = -i{df,9iA - aif, f - i{d^,6iB + ai^^f 

-iid^02A - a2M)' - i{d^02B + a2^,f 

+gcos{eiA + 0iB+O2A + d2B)- (38) 

If there is no gauge field, the condensation of 9s would 
lead to four gapless Goldstone modes. However, in the 
condensate, mode 0ia + Oib + 02A + d2B is gapped out 
by the g term in (j38p . this implies that in the superfluid 
phase 6iA + ^is = -(^2A + ^2b)- Meanwhile, Oia ~ Oib 
will gap out flip through the Higgs mechanism, 62A — 
62B will also gap out through the Higgs mechanism, 
therefore the only gapless mode in the condensate is 9ia+ 
QiB = —{O2A + Q2b)- 

Notice that S} can create a pair of V'lA and ipiB par- 
ticles and also can annihilate a pair of 'ip2A and "025 par- 
ticles. Therefore we can identify 

^ eMKOiA + Oib)] - exp[-z(02A + e2B)]- (39) 

Thus, the Goldstone mode Oia + Oib is exactly the global 
U{\) phason mode of ~ exp(z6'). 

If one approaches the transition from the superfluid 
phase, the transition can be viewed as condensation of 
vortices in the superfluid. There are four components 
of vortices, corresponding to the four flavors of matter 
fields. The gapless Goldstone mode becomes the non- 
compact C/(l) gauge field in the dual language. The ver- 
tex operators in the height field language are the vortex 



tunnelling terms. The vertex operators create or annihi- 
late gauge flux of the original gauge fields oi^ and 02^. 
For instance, exp(27rz(pi) creates one unit flux of ai^. As 
pointed out in references^SiSZ,, when one flavor of gauge 
field is coupled to two different matter fields, the vor- 
tex of each matter field carries half flux quantum. Since 
vortex Via and vib carry opposite gauge flux, the ver- 
tex operator cos(67r(y9i) corresponds to tunnelling process 

v^^a'^Ab + H.c The dual Lagrangian can be effectively 
written as 

L - -t|(a^ - iA^,)viA? - t|(a^ - «A^)fiB|' 
-t|(a^ + iA^)v2A? - t|(a^ + «A^)f2s|' 

+9{vfAvlB + H.C.) + g{vl\vlB + H.c.) 

+g2{v\j^ViBvlgV2A + H.c.) + ri{viAViBV2BV2A + H.c.) 

+5i(w1^wisw|^w|b + H.c.) + gi{vlAV2Bvf'AvlB + H.c.) 

(40) 

in is the dual form of the Goldstone "phason" 
mode in the superfluid phase, g, gi and 32 are the tun- 
nelling terms due to the vertex operators in (j23[) . Tun- 
nelling term rj is independent of monopoles, as this term 
conserves the total vorticity (consistent with the U{1) 
gauge symmetry of (|40p). and also conserves the total 
gauge flux of the gauge fields ai^ and a2fi, therefore it 
should exist in the field theory. 

Let us denote Via ~ exjp{—ixia)- After the condensa- 
tion of vortices, modes xiA — XiB and X2A — X2B are 
gapped out by the monopoles. Mode xiA + XiB + X2A + 
X2B are gapped out by the 7 term in equation ([3D]), i.e. 
XiA + XiB = -{X2A + X2b); also xiA + XiB is gapped 
out by A^ through the Higgs mechanism. Therefore in 
the condensate of vortices, there is no gapless excitations, 
which is consistent with the crystalline phase. 

If the monopole effect is turned off, the transition point 
is described by two gapless noncompact gauge fields and 
four flavors of matter fields. However, whether these gap- 
less gauge fields and matter fields can survive when the 
monopoles are turned on is an open question. If the 
monopoles gap out the gauge field and confine the mat- 
ter fields, at the transition there is no gapless excitation. 
In this case our theory predicts a direct first order tran- 
sition. 

The superfluid phase and the plaquette phase break 
different symmetries, and according to the classic Landau 
phase transition theory, the transition should be either 
flrst order, or split into two transitions, with a disordered 
phase (or a phase with both orders) in between. There 
is no universal law that guarantees one direct first order 
transition. 

In our theory, the intermediate phases can be under- 
stood as the condensate of composites of defect ij^ia. 
A gapped disordered phase can be obtained if compos- 
ites which only carry local gauge charges but no global 
U{1) charge are condensed. For instance, if composites 
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Wia'^ib and iplj^ip2B are condensed while all the other 
composites are disordered, the gauge fields are gapped 
through the Higgs mechanism, therefore the height fields 
are disordered, the crystalline order disappears. Also, 
since the composites carry zero global U{1) charge, there 
is no gapless Goldstone mode. Thus we can conclude 
that the condensate of V'iaV'is and •4'2a'^^b ^ spin 
disordered phase i.e. a spin liquid phase. Notice that 
iplj^ipiB ("02^^25) carries two unit gauge charges of 
gauge field ai^ {0'2^l), therefore the condensate of ip\y^tpiB 
and ipl^ip2B is a spin liquid with Z2 x Z2 gauge sym- 
metry, which is the residual gauge symmetry after the 
condensation of the composites of matter fields. On 
the other hand, if composites carrying only global U{1) 
charge are condensed, the superfluid order should coex- 
ist with the crystalline order. For instance, the compos- 
ite '!/'iaV'ibV'2aV'2B does not carry any gauge charge, the 
condensate of this composite is a superfiuid order, and 
the crystalline order still exists. Thus this phase is a 
supersolid phase. 



VI. TRANSITION TO THE NEMATIC PHASE 

The existence of nematic phase can be derived easily at 
the negative large D limit. When D is negative and be- 
comes the dominant term in the Hamiltonian ([T]) , the sys- 
tem is effectively a spin- 1/2 system, since on each site 
can only be ±1. The classical ground state of this model 
is that every unit triangle should have either (1,1,-1) or 
(1,-1,-1) configuration. This classical ground state is 
the same as the classical Ising model on the kagome lat- 
tice, with large degeneracy. If the same Boltzman weight 
is imposed for each classical ground state, the kagome lat- 
tice Ising model is disordered, and the correlation length 
is finite^S-. By contrast, a related classical system is the 
classical Ising model on the triangular lattice, while if 
the same Boltzman weight is imposed for each classical 
ground state, the triangular lattice Ising model is critical, 
and an infinitesimal quantum perturbation is relevant at 
this critical point and drive the system into a crystalline 
phase. If infinitesimal transverse magnetic field ha^ is 
turned on, the system is driven to V3 x Vs order^^; if 
ferromagnetic exchange — Jxy{(jf aj +af (Jj) is turned 
on, the system is driven into a supersolid phase, which 
breaks both U{1) symmetry ( (cr+) 7^ ), and transla- 
tional symmetr}*2^. Unlike the Ising model on the tri- 
angular lattice, the classical Ising model on the kagome 
lattice is disordered, with finite correlation length. In the 
original Hamiltonian ([1]), if \Jxy\ ^ D, the second order 
perturbation generates a term which fiips — 1 state to 
5^ = — 1 state and vice versa. The effective Hamiltonian 
reads 

Heff = E + -^M) + J^^^^r (41) 

<»,i> 



t ^ J^y/D. As studied inSi2£^ for the spin-1/2 system on 
the kagome lattice, infinitesimal ferromagnetic XY ex- 
change yields superfluid order, (o-^) 7^ 0. The spin-1/2 
raising operator is the nematic order parameter (5^)^, 
therefore infinitesimal \Jxy\ drives the system into a ne- 
matic phase. 

Although the nematic phase and the plaquette phase 
do not necessarily touch each other in the phase diagram, 
a direct transition between these two phases is possible 
when they are adjacent in the phase diagram. It is con- 
ceivable that a certain type of spin Hamiltonian can real- 
ize the direct transition between the nematic phase and 
the six-fold plaquette phase. This direct transition is 
more likely to occur when Jxy > than the case with 
Jxy < 0. In the case with Jxy > 0, every hexagon is 
effectively penetrated by one tt fiux of oi^ and one tt flux 
of a2fi- The motion of defects is strongly affected by the 
background magnetic fields, and several interesting pos- 
sibilities can happen. One of the possibilities is that the 
defects condense in pairs, i.e. ((i/'ia)^) 7^ 0, as a pair 
of defects does not see any background flux. After the 
pair condensation, the Goldstone mode is 2{9ia + Oib), 
corresponding to the phase of (S"^)^, so the system is 
in the nematic phase discussed above^. One important 
difference between the nematic phase and the superfluid 
phase is that, each vortex in the nematic phase only car- 
ries one quarter fiux of the gauge fields, therefore the 
vertex terms in (j23p correspond to even higher order of 
vortex tunnelling processes. 

A direct transition between the nematic phase and the 
plaquette phase can be described by the following action 
of paired matter field ^'ia = (V'ia)^ 

L = -t\{d^ ~ 2ia,^)^iA\^ - t\{d^ + 2mi^)*iBp 

-t\{d^ - 2ia2^)^2A\^ - t\id^ + 2m2M)*2B|' 

-Hi + (42) 

Again the ellipses include the monopole terms, and Hj 
contains all the possible interaction terms between mat- 
ter fields. Just like the four-defect creation term dis- 
cussed in the previous section, 5^'iyi^is^2yi^2S + H.c. 
should in principle exist in the interaction, thus phason 
mode X]?=i ^a=A gapped out in the condensate of 
"^ia- Without the monopole terms this transition is a 
gapless second order transition. 

Now the question boils down to if the monopole effect 
is going to be relevant at the critical point described by 
action (|42p . Since the nematic phase is a pair conden- 
sate, each single vortex in the nematic phase carries only 
one quarter fiux of each flavor of gauge fields, so the ver- 
tex operator in (j23p corresponds to even higher order of 
tunnelling processes than the superfiuid case. The dual 
action now reads 
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-t\{d^ + iA^)v2A\^ - t\{df, + a^)w2s|' 
+9i4>tB + H.c.) + givl'A + H.c.) 

+ H.c.) + r]{viAVlBV2BV2A + H.c.) 

5l(^'^A^lB^'2A«2B + H.c.) + giivl^vlgvWvfg + H.c.) 

(43) 

g, gi and 172 terms are vortex tunnelling processes cor- 
responding to the vertex operators in ()23p . notice that 
now Via and vib iv2A and V2b) carry one quarter unit 
flux of icL2fj.) ■ V term is a tunnelling which does not 
rely on monopole, as it not only complies with the U{1) 
gauge symmetry of the dual action but also con- 

serves the flux numbers of the original gauge fields oi^ 
and 02/^. Following the similar argument as referencesii^, 
the monopole terms (vortex tunnelling terms) are likely 
(but not rigorously proved) irrelevant at the transition 
fixed point. It is known that at the 3D XY fixed point, 
Z4 anisotropy is irrelevant, while the Z3 anisotropy as in 
could be relevant. However in our case, the vertex 
operators in (j23p could be irrelevant at the transition due 
to the pairing of gauge charges ((^ . Although the vertex 
operators are always relevant at the Gaussian fixed point 
of (IM)) , it could be irrelevant at the order-disorder tran- 
sition of height fields. It is expected, as in the calculation 
that follows, that the scaling dimension of the vertex op- 
erators is approximately proportional to the number of 
flavors of matter flelds, and proportional to the square of 
the product of electric charge and magnetic charge'^'^. 

We can roughly estimate the scaling dimension of the 
monopole operators from a random phase approximation 
(RPA) calculation. After integrating out the Gaussian 
part of the matter flelds in ([^^ . an effective action for 
gauge flelds ai^ and a2^ is generated 

2 . ,3, 

^=E J j^Nn'\k\<j,\a{k)o.\' + (44) 

N is the number of flavors of bosons coupled to each 
gauge fleld, n is the number of gauge charge carried by 
each boson. In our case N — n — 2. In the dual theory, 
the kinetic term for the height fleld is softened to be 
~ fc^, and the monopole energy diverges logarithmically 
instead of converging in the infrared limit^i^-^ . The dual 
height flelds now have the action 

a— 1 

From this calculation one can see that the scaling dimen- 
sion of vertex operators is proportional to Nn'^ . 

Hy in contains three types of terms. The scaling 
dimensions for coslGiripa) and cos(27r(y9i -I- 47r(p2) calcu- 
lated from the RPA approximation is higher than the 
Z4 anisotropy studied before^'^. The third vertex opera- 
tor is cos(27r((^i — 1^92)), the scaling dimension calculated 



from RPA is higher than the Z3 anisotropy of 3D XY 
flxed point, also, on the RPA level, the scaling dimen- 
sion is equal to the case with Z4 anisotropy and iV = 1 
discussed in reference^, which has been shown to be irrel- 
evant at the transition between the Higgs phase and the 
conflned phase. Recently a Monte Carlo simulation has 
shown that the transition between the crystalline phase 
and the superfluid phase in a bosonic model with 1/3 
filling on the kagome lattice is a very weak first order 
transitions^, on the RPA level, the scaling dimension of 
the monopole in that case is smaller than the dimension 
of all the triple vertices and very close to the scaling di- 
mension of cos(27r((/3i — 1^2)) in our case. Therefore it is 
possible that the vertex operators in our problem are ir- 
relevant at the transition between the nematic phase and 
the plaquette order. When the vertex operators are irrel- 
evant, the critical point is a direct gapless second order 
transition, with four flavors of deconflned matter fields, 
as well as two flavors of noncompact gauge fields. 

VII. TRANSITION TO THE ^3x^3 PHASE 

When J„j > and much larger than other coefficients, 
the state is most likely to be either the \/3 x \/3 order 
in Fig. [U or the g = state in Fig. H From the 1/S 
expansion, this VS x \/3 order is supposed to be the 
global minimum of all the classical degenerate ground 
states of the Heisenberg model on the kagome lattice at 
the isotropic point^^, although the q — state has also 
been proved to be one of the local minima. Both \/3 x \/3 
and q — states are very typical configurations for spins 
on kagome, they can be stabilized by the second or third 
nearest neighbor interactions. Also, since both states are 
coplanar, they are expected to be even better candidates 
in the large Jxy case. 

Since now the defect hopping is frustrated by the back- 
ground magnetic flux of gauge field ai^ and 02^ through 
each hexagon, the phase angle of the defects cannot be 
uniformly distributed on the whole lattice. We will see 
that the \/3 x \/3 phase can be interpreted as the conden- 
sate of the four flavors of charge fields in the background 
gauge fluxes. 

Because of the interaction between different matter 
flelds g{ipiA4'iB^2Aip2B + H.c), we have the following 
relation between the phase angles 

0^dlA + OiB = -{O2A + d2B). (46) 

6 is the phase angle of . The distribution of phase 
6 can be deduced from the distribution of 9ia and 9ib- 
Notice that tpiA and tpis both live on the sites of the 
honeycomb lattice, and hop on two different triangular 
sublattices (Fig. [8]). With a background magnetic field 
ai^, the effective Hamiltonian for the motion of ^ia is 

H= J2 tcos{0iA.~OiA^j) + (47) 

<i,j> 
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FIG. 8: Defects which violate the 3-color constraint hop on 
one of the two triangular sublattices (red and green) of the 
honeycomb lattice. 




FIG. 9: The \/3 x ^/3 order of 9 obtained from the ordered 
pattern of 9ia and 9ib- The black arrow corresponds to angle 
9, and the pattern of 9 can be obtained from its adjacent red 
arrow (6*1,4) and green arrow {9\b)- 

This is an antiferromagnetic XY model on the triangular 
lattice, and after the condensation of Oia^ the ground 
state is the VS x \/3 order. This phase can be viewed 
as the staggered vortex density phase on the triangular 
lattice. 

If both tpiA and ipiB condense (due to the time rever- 
sal symmetry, if ipiA and ipiB condense, 'ip2A and "023 
will also condense), the phase angle 9 can be determined 
from the distribution of 9xa and 9ib- By adding the two 
ordered patterns of both 9ia and Oib together, the or- 
dered pattern for 9 is automatically obtained, and the 
order can only be either g = state or the -s/S x ^/2> state 
(Fig. and (Fig. fTU]) . In these ordered phases, the 
Goldstone mode is still 9 — 9ia + 9ib- 



VIII. LONGITUDINAL MAGNETIC FIELD 

A longitudinal magnetic field along z breaks the time 
reversal symmetry in the Hamiltonian, and much of the 
physics is significantly changed. Let us assume the mag- 
netic field is small, i.e. /iS*^ is much smaller than 




FIG. 10: q = Q order of 9 obtained from ordered pattern of 
9iA and 9\b. 

and D in the model. Note that this precludes accessing 
strong-field phenomena such as magnetization plateaus 
in this theory. 

The six-fold degenerate plaquette phase is expected to 
survive in a small longitudinal magnetic field. Since in a 
small magnetic field, the classical ground states without 
Jxy are still configurations with (1,-1,0) triangles only, 
and therefore the ring exchange term generated by J^y 
is still going to select the plaquette ordered state as the 
ground state. 

However, the excitation energies of defects are changed 
by the longitudinal field: (1,1,-1) defects have lower 
gap than the (1,-1,-1) defects. Therefore when |Ja:j,| 
is turned on, (1,1,-1) defects should condense before 
(1,-1,-1). As we will see in this section, the conden- 
sate of (1,1,-1) defect is actually a supersolid phase, 
with both global U{\) symmetry breaking and the space 
symmetry breaking. 

Let us take the case with J^y < as an example. As 
long as (1,-1,-1) defects remain confined and gapped, 
the total number of (1,1,-1) defects is conserved be- 
fore its condensing, because (1, 1,-1) defects cannot be 
excited without (1,-1,-1) defects due to the conser- 
vation of total 5^. Therefore when (1,1,-1) defects 
condense and (1,-1,-1) defects remain confined and 
gapped, the system still has a gapless Goldstone mode 
due to the spontaneous breaking of the global conserva- 
tion of (1,1,-1) defects. The gapless Goldstone mode 
manifests the superfluid phase. 

Secondly, the spatial symmetry breaking can still be 
studied in terms of the dual height fields. Because 
(1, 1, —1) is the vortex of height field the condensate 
of (1,1,-1) charge is the disordered phase of (pi, ipi no 
longer has nonzero expectation values. However, because 
Lp2 is still ordered, the vertex operator —acos{(!)'KLp2) 
in (j23p has 3 minima, corresponding to 3 fold degener- 
ate states. These 3 minima are the plaquette orders of 
(—1, 0, —1, 0, —1, 0) hexagon on 3 different sublattices. 

The height fields Lpi and Lp2 are coupled through the 
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vertex operators as shown in (j23p . This coupHng is not 
going to Hft the 3 fold degeneracy of (^2 when ipi is dis- 
ordered. The reason is as follows: The whole action 
(|24p is invariant under transformation (pi —^ ipi + 1/3, 
V2 ^ </32 + 1/3. Since (pi is disordered, after integrating 
over ifi, the effective action Heff{(p2) for ip2 does not 
break this symmetry, and the leading vertex term gener- 
ated from integrating out ipi is cos(67r(p2)- This can be 
clearly seen from the following equations 

exp(-i?e//(<^2)) = J D(fi exp(-iJ(.^i, (^2)) 

= J Dipt exp(-iJ(^i + 1/3, (^2 + 1/3)) = 

j + 1/3) cxp(-i/(^i + 1/3, (^2 + 1/3)) 

= j i?(^iexp(-iJ(^i,(^2 + l/3)) 

= exp(-i/e//(<^2 + l/3)). (48) 

Notice that the above proof is only valid if Lpi does not 
take any nonzero expectation value, i.e. (pi is in disor- 
dered phase. Alternatively, one can understand this ar- 
gument from the conservation of the gauge fluxes. Vertex 
operator cos(27riy9Q,) can annihilate or create one unit flux 
of gauge field Cq^. The total flux of both gauge fields is 
conserved mod 3 in vertex operator Hamiltonian (j23p . 
and as the disordered phase of tp^ (the condensate of 
(1, 1,-1) defect) does not tend to violate this conserva- 
tion, the resultant effective Hamiltonian after integrating 
out ipi fields does not break the Z3 conservation of total 
gauge flux, i.e. the lowest order vertex operator of the re- 
sultant effective Hamiltonian of Lp2 is cos(67r(/32)- There- 
fore the 3-fold degenerate plaquette order is not lifted. A 
similar result is obtained for J^y > too. Thus, the phase 
with defect (1, 1, —1) condensed while defect (1,-1,-1) 
confined breaks both spatial symmetry and the global 
U{1) symmetry, and hence must be the supersolid phase. 

When D <Q and large, a small longitudinal magnetic 
field changes the physics severely. In this regime, the 
classical ground state has Ising configuration (1,1,-1) 
on each triangle. In the previous sections we mentioned 
that the classical Ising ground state on the kagome lat- 
tice is disordered with finite correlation length. However, 
once the longitudinal magnetic field is turned on, every 
triangle has (1,1,-1) configuration. Since the sites of 
the kagome lattice are the links of the dual honeycomb 
lattice, the ground state configurations with small mag- 
netic field can be mapped onto the dimer configurations 
on the honeycomb lattice, with = —\ mapped onto 
dimer, and 5^ = 1 mapped onto empty link. If the same 
Boltzmann weight is imposed on every dimer configura- 
tion, the system is again critical, with power law decaying 
spin-spin correlation function. 

Since the classical ground state is critical, it is again 
very instable against quantum perturbations. If a small 
Jxy is turned on, the system is driven into a gapped crys- 
talline phase. At the sixth-order perturbation, a ring 




FIG. 11: The distribution of 5^ = 1 and 5"^ = -1 sites on 
the kagome lattice is equivalent to the distribution of dimers 
of dimer model on the dual honeycomb lattice. The ring ex- 
change term (I49p plays the same role as the dimer resonating 
term in the quantum dimer model. 

exchange term is generated 

H„ng = J2 ~tiiSlS2)HslSi)^iSlSef + H.c.) (49) 
o 

t ~ J^y/ [ JID'^). This ring exchange term plays the same 
role as the dimer flipping term in the honeycomb lat- 
tice quantum dimer model, which will generally lead to 
a crystalline phase. Notice that, besides the off-diagonal 
flipping term in (|49p . diagonal terms are also generated. 
According to several previous works^'''^^, the diagonal 
terms generated by perturbation theory favor flippable 
hexagons, therefore it is expected that the crystalline or- 
der is either plaquette order or columnar order-?^. 

Presumably the crystalline phase disappears when 
J^y/D ^ h, Since the sign of t in (1^5]) is always positive 
(independent of the sign of Jxy), the crystalline phase 
should extends symmetrically on the two sides of the 
classical line with J^y = 0, until the system enters the 
nematic phase. The sketchy phase diagram in a small 
magnetic field is shown in Fig. [121 note that this phase 
diagram involves a lot of phases, the detailed topology 
of the phase diagram would depend on the details of the 
microscopic model. 



IX. THE SU{2) POINT 

At the isotropic point Jxy — Jz and D — 0, the v^x \/3 
state is just one possibility. This state is obtained from 
quantum perturbation on the classical limit. If we start 
with the quantum limit, another possible state can be 
obtained: the q = dimer plaquette state. 

This state can be understood quite easily from the 
quantum dimer model on the kagome lattice. For a spin- 
1 system, each site can form two spin singlets, which 
means each site connects to two dimers. Since every site 
on the kagome lattice is shared by four links, this dimer 
model is half filled. The dimer resonance term on the 
kagome lattice is shown in Fig. [121 which can flip the 
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FIG. 12: The phase diagram in smaU longitudinal magnetic 
field h. The difference between this phase diagram and the 
one without magnetic field in Fig. [5] is that, between the 
plaquette phase and the superfiuid phase there is a supersolid 
phase. And in the negative and large D case, small J^y is 
going to generate a gapped crystalline order. 




FIG. 13: The dimer model on the kagome lattice. Since this is 
a spin-1 system, every site connects to two dimers. The dimer 
hipping term flips configuration Ci to C2 and vice versa. 



dimer covering Ci to C2 and vice versa. The dimer model 
Hamiltonian is now written as 

H = -i(|Ci)(C2| + H.c.) + V{\C^){C^ + |C2)(C2|). (50) 

As long asV <t, the exact ground state wave function 
of this Hamiltonian can be written as 



I*) 



o 



IC2)). 



(51) 



This state does not break any space symmetry, and it 
minimizes the energy of each hexagon individually. This 
state should be the ground state of a certain type SU (2) 
invariant spin Hamiltonian. 

Now the question is whether this dimer plaquette phase 
is a new phase or it can be continuously connected with 
one of the other states discussed early this paper without 
any physical singularity. Notice that the gapped photon 
phase with (0, 0, 0) on every triangle breaks no space sym- 
metry either, thus one can imagine adding D{S^')^ on the 



V 



FIG. 14: A typical state in the Haldane phase. In the bulk, 
every site is shared by two dimers. However, at each edge, 
there is one residual unpaired spin-1/2 variable. 




FIG. 15: The edge of AKLT state on the square lattice. There 
is one extra spin-1/2 degree of freedom on each site of the 
edge, therefore the edge state is effectively a spin-1/2 spin 
chain, which is either gapless or breaks translational symme- 
try. The AKLT state is qualitatively different from the state 
with = everywhere. 



isotropic Hamiltonian, and the ground state wave func- 
tion can be continuously deformed to the gapped photon 
phase. 

One has to be careful about the naive argument above. 
Let us consider the one dimensional analogues of the 
dimer plaquette phase and the (0,0,0) phase as a check 
of our naive argument. The one dimensional Haldane 
phase^I for spin-1 Heisenberg chain is gapped, and breaks 
no symmetry. One can imagine that by adding D(S'^)^ 
on the Heisenberg Hamiltonian the Haldane phase will be 
continuously connected to the state with 5*^ = every- 
where. However, the spin-1 Heisenberg chain is charac- 
terized by two special properties: the first is the existence 
of gapless edge states, the second is the hidden diluted an- 
tiferromagnetic order. The existence of the gapless edge 
states can be understood as follows: All the sites in the 
bulk are shared by two dimers, while the site at the edge 
only connects to one dimer, hence there is a residual spin- 
1/2 degree of freedom at each edge (Fig. The hidden 
diluted antiferromagnetic order can be viewed from ex- 
panding the AKLT state (the explicit wavefunction of 




FIG. 16: The edge of the kagome lattice. Every unit cell at 
the edge has even number of spin-1/2 quantities. Presumably 
the edge state is gapped and featureless. 
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one particular state in the Haldane phase) of the spin-1 
chain in the basis of . In this expansion, one typical 
state is as follows 

+ -00000 + 000 -+00 - +000 (52) 

Every = 1 site is always followed by one ~ —1 
site, although there could be a number of 5*^ = sites 
in between. A special nonlocal string operator could be 
introduced to describe this hidden order in the Haldane 
phase^^i^. 

The nice features of the Haldane phase also exist in 
the AKLT state of spin-2 systems on the square lattice. 
Let us take a cylinder geometry as an example. There 
is one unpaired spin-1/2 degree of freedom on each site 
of the edge, therefore the edge state is effectively a spin- 
1/2 chain, which is either gapless or gapped but breaks 
translational symmetry. The AKLT state is qualitatively 
different from the state with S*^ = everywhere as well 
(Fig. [HI). 

However, the edge states are missing in the dimer pla- 
quette phase of the spin-1 dimer model on the kagome 
lattice. If we take a cylinder with edges, on each site 
of the edge there is also a residual spin-1/2 degree of 
freedom. However, due to the geometry of the kagome 
lattice, in every unit cell of the edge there are even num- 
ber of spins, therefore effectively the edge state is a chain 
with integer spin (Fig. [TB|) . The resultant edge state is 
generally gapped and featureless at the edge, which is 
the Haldane phase on a closed circle. Thus, we conclude 
that the dimer plaquette phase can be continuously con- 
nected with the gapped photon phase, with (0,0,0) spin 
configuration on every triangle. 

X. OTHER TRANSITIONS 

In the phase diagrams Fig. [5] and Fig. [T^l there are 
several other transitions which are interesting. First of 
all, in Fig. O the transition between the nematic phase 
and the supcrfluid phase is probably an Ising transition, 
as this transition breaks the I2 symmetry in the nematic 
state. The transition between the nematic phase and the 
•\/3 X \/3 phase is supposed to be a first order transition. 

In the case with magnetic field, since another crys- 
talline order is opened up (Fig. I12p . there is a transition 
between the crystalline phase and the nematic phase. 
However, now that in the case of large and negative 
D, the system can be described by an effective spin-1/2 
model (|4T|) , the transition can be understood as the tran- 
sition between the crystalline order and the superfluid 
order of hard core bosons on the kagome lattice. This 
transition has been studied in reference a^'^'^" , and the 
transition is a weak first order transition. 



XI. CONCLUSIONS AND EXTENSIONS 

In the current work we studied the global phase dia- 
gram of the spin-1 XXZ antiferromagnet on the kagome 
lattice. Various phases which have been studied before 
can be obtained from condensation of the defects in one 
single phase. The phase diagram was also obtained for 
the case of a magnetic field along the z direction. One 
route to test this phase diagram experimentally is by 
neutron scattering or other measurements on the spin- 
1 kagome materials, for instance Ni3V208 and BF4 salts. 

In all the previous sections, the model under consider- 
ation only contains quadratic interactions. However, in 
some circumstances, for instance a spin-1 bosonic system 
trapped in an optical lattice, the biquadratic interactions 
— J2{Si ■ S2Y have been shown to be important. This 
biquadratic term can help to stabilize the nematic phase, 
when it becomes the dominant term in the Hamiltonian. 
In closing we briefly explain one interesting consequence 
of this biquadratic interaction, in case a cold-atom real- 
ization of this Hamiltonian is constructed. 

Suppose that the system is in the six-fold degenerate 
plaquette ordered state, and let us gradually turn on the 
biquadratic term in the XY plane 

Hu = -J2{S^S^ + SySjf. (53) 

This biquadratic term is consistent with the XXZ sym- 
metry of our model. At the third order in pertur- 
bation theory, this biquadratic term can generate res- 
onance between (+1, — 1, +1, — 1, +1, — 1) hexagon with 
(— 1, +1, — 1, +1, — 1, +1). Notice that, although spin 
— +1, —1 and are treated as three colors, the Z3 
symmetry is missing in Hamiltonian ([9]), since the reso- 
nances were only between (+1, 0, +1, 0, +1, 0) hexagons 
and between the (— 1, 0, — 1, 0, — 1, 0) hexagons. There- 
fore the full Z3 symmetry can be restored by turning 
on the biquadratic term ([53)) . At this Z3 point the 
ground state is probably a nine-fold degenerate plaque- 
tte order. Once the biquadratic XY exchange dominates 
the quadratic XY exchange, the phase becomes a three- 
fold degenerate plaquette ordered state with resonating 
(+1, —1) hexagons. 
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